Onset of superconductivity and hysteresis in magnetic field for a long cylinder 
obtained from a self-consistent solutions of the Ginzburg— Landau equations 
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Based on the self-consistent solution of a nonlinear system of one-dimensional GL-equations, the 
onset and destruction of superconductivity, the phase transitions and hysteresis phenomena are 
discussed for a cylinder (radius R) in an axial magnetic field (H) for arbitrary R, X, H, m (x is 
the GL-parameter, m is the total vorticity of the system). The edge-suppressed solutions (which 
are connected with the jumps of magnetization in the states with fixed vorticity m), the depressed 
solutions (responsible for the hysteresis in type-II superconductors), and the precursor solutions 
(which describe the onset of superconductivity in type-I superconductors) are also studied. The 
limits of applicability of the so-called linear equation approximation are discussed. 
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I. INTRODUCTION 

The macroscopic GL-theory [1] is widely used for 
studying the behavior of superconductors in magnetic 
field. This theory leads to two coupled nonlinear three- 
dimensional equations for a complex order parameter \t 
and magnetic field vector-potential A. Because, in a 
general case, it is impossible to find the exact solution 
of this system analytically, the numerical computations 
are used, or various limiting particular cases are stud- 
ied. Thus, the London theory [2] follows from the GL- 
equations [1] in the limit ip = \*5f\ — 1, and the Ginzburg 
approach [3] corresponds to the case of a constant or- 
der parameter, which depends on the weak external field 
H, ip = tp(H). Another approximation was formulated 
by Saint- James and de Gennes [4], who assumed, that the 
nucleation of superconductivity in large fields is described 
by small values of the order parameter (ip = \^f\ <C 1). 
Instead of full system of nonlinear GL-equations, they 
considered a single linear one-dimensional equation for 
ip, what allowed them to find analytically the nucleation 
field H c3 = l.69H c2 (H c2 = (p / '(2tt£ 2 ), £ is the coher- 
ence length) for a bulk superconducting slab in paral- 
lel magnetic field. Later, Saint-James [5] used the same 
approach to find the oscillation dependence of the nu- 
cleation field H C 3 versus applied field H for a very long 
circular cylinder of finite radius R. In the papers [4,5] 
the basic conception of the surface superconductivity was 
laid down, according to which the onset of the supercon- 
ductivity begins at the fields H ~ H c3 in a form of thin 
superconducting sheath, layered near the surface of the 
specimen (the so-called "giant-vortex state"). 

The structure of the giant-vortex solutions with finite 
values of the order parameter (ip ~ 1) was studied nu- 
merically in a number of scattered publications (for the 
review of early works see [6]). A more systematic study 
of such solutions, carried out recently in [7], has revealed, 



that the axially symmetric solution of fixed vorticity m 
has different form, depending on the value of the in- 
creasing external magnetic field H (m = 0, 1,2- ■■ is 
the quantum number, ensuring that \& is single valued). 
In type-II superconductors the giant-vortex solution has 
the usual Meissner-type space-profile for fields H < Hi 
(where Hi(m, R, x) is some critical field). At H = Hi 
this form becomes unstable, and for H > H\ the giant- 
vortex solution acquires (by a first-order jump) a new 
"edge-suppressed" form [7]. With the field H further in- 
creasing, the edge-suppressed solution degenerates grad- 
ually and vanishes finally by a second-order phase tran- 
sition (ip — > 0) at some field H2(m,R,x). (The field 
H2(m, R, x) is just the critical field, found in [5] from a 
linear theory; for R ^> £ and m ~ 1 the field H^ coincides 
with H C 2] for R 3> £ and m 3> 1 the maximum value of 
the field H 2 coincides with H c3 — 1.69H c2 .) 

In type-I superconductors the destruction of the super- 
conducting state occurs at the field Hi(m, R, x) (i.e. at 
the field of maximal superheating) by a first-order jump 
from s-state (having a finite value of ip ~ 1) to n-state 
(ip = 0). Because in this case there are no solutions 
with ip 1, the linear theory is not applicable, and the 
superheated s-boundary, found from the full system of 
equations, deviates substantially from that, predicted in 
[5]- 

As is shown in the present paper, the nucleation of 
superconductivity from n-state in decreasing fields also 
proceeds differently in type-I and type-II superconduct- 
ing cylinders. 

In type-II cylinders the nucleation begins at 
H2(m, R, x) as a second-order phase transition from n- 
state, then a small nucleated s-solution grows, repeating 
the edge-suppressed form down to the field H\(m, R, x). 
At this field the shape of the order parameter starts to de- 
viate gradually from the edge-suppressed form and passes 
continuously into a new "depressed" form, which exists 
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down to some "restoration" field H r (m, R, x), where the 
first-order jump to the Meissner-type solution occurs. In 
the field interval A = H \ — H r there exist simultane- 
ously two solutions for the same field H: one is the usual 
Meissner-type solution (ip s=a 1), and the second is the de- 
pressed solution (ip < 1), which describes the hysteretic 
s-state of the cylinder. The region of parameters, where 
A(m, R, x) > (i.e. the phase boundarywhere the hys- 
teretic transitions between superconducting states of the 
same vorticity may exist), is found self-consistently. 

In type-I cylinders (in the field decreasing regime) a su- 
percooled n-state persist down to some field H p , where 
the feeble (-0 <§C 1) "precursor" solution forms, what cor- 
responds to the second-order (n, s)-phase transition, so 
the field H p may be found from linear equation approx- 
imation [5]. However, the precursor solution exists only 
in a very small interval (A p ) of fields below H p (the field 
interval A p = H p - H r ~ 1CT 2 - HTfHp), after that 
the first-order restoration of the full Meissner solution 
(ip ss 1) occurs. [Thus, the nuclcation of superconduc- 
tivity in type-I cylinders is an "almost first-order" phase 
transition.] The phase region for type-I cylinders, where 
the magnetically " supercooled" normal state [3] may ex- 
ist, is also found self-consistently. 

As was mentioned above, there are two different 
regimes to study the effect of the external field penetra- 
tion into the superconductor interior: the field increase 
and the field decrease regimes, which should be consid- 
ered separately. It is demonstrated below, that if the field 
increases, there exists a critical line S\-u on the plane of 
variables (to, R, x), which separates two superconducting 
regions, s\ and s\\. In the region s\ a superconducting 
state vanishes by a first-order phase transition to the nor- 
mal state (by a jump from ip ^ 0). The superconductors, 
which belong to the si-region, is naturally classified as 
type-I. (This classification is applied to one-dimensional 
centrally symmetric states of a cylinder.) In the region 
sn (if the field increases) the destruction of superconduc- 
tivity is a second-order phase transition to normal state 
(ip — > 0). The superconductors, which belong to the sn- 
region, naturally classified as type-II. 

It is shown below (see also [7] ) , that the critical bound- 
ary 5*i_ii between type-I and type-II superconducting 
cylinders is not x c = 1/ \/2 (as in the case of normal and 
superconducting half-spaces, brought into contact [8]), 
but is a complicated function of the cylinder radius R 
and vorticity to, i.e. x c = f(R,m). It is shown, that for 
sufficiently small R, all superconducting cylinders (inde- 
pendent of the GL-parameter x) belong to type-II class. 

The paper is organized as follows. In Sec. II the basic 
GL-equations are written, as well as the linear equation 
for tp. The known analytic solutions of the linear equa- 
tion are reproduced (in terms of the Kummer functions). 
The nonlinear generalization of the Kummer-type equa- 
tion is regarded. In Sec. Ill the phase diagrams, which 
divide (on the plane of variables (R, x)) the regions of 



first- and second-order phase transitions to the normal 
state (and also the regions, where the magnetic hystere- 
sis is possible), are presented. The difference between the 
depressed and edge-suppressed solutions is explained. 

In Sec. IV the critical fields, found self-consistently 
from the GL-equations for a cylinder in the increasing 
field, are compared with those, found from the linear ap- 
proach, for both type-I and type-II superconductors. In 
Sec. V the analogous comparison is made for the de- 
creasing field. The "supercooling" of the normal state, 
the precursor solutions and hysteresis in type-I supercon- 
ducting cylinders are also discussed. Sec. VI contains a 
short resume and discussion of the results with possible 
connection to experiment. 



II. EQUATIONS 

In what follows below, the case is considered of a long 
superconducting cylinder of radius R, in the external 
magnetic field H > 0, which is parallel to the cylinder 
element. Only the radially symmetric (one-dimensional) 
solutions are studied. In the cylindrical co-ordinates the 
system of GL-equations may be written in dimensionlcss 
form [7] 
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The dimensioned potential A, field B and current j s are 
related to the corresponding dimensionless quantities by 
the formulae: 
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The total vorticity to in (3) specifies how many flux 
quanta are associated with the vortex, centered at the 
cylinder axis (the so-called giant- vortex state [5,9]). [The 
equivalent names for m, are: vorticity, fluxoid, orbital 
momentum, magnetic quantum number, winding quan- 
tum number.] 

The boundary conditions to Eq. (1) are: 
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where i? e = h 6 = H/H e , H € = o /(2< 2 ). 
The boundary conditions to Eq. (2) are: 
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(We use here the coherence length £ as a unit of length, 
instead of the field penetration length X — x^, used in [7]. 
Both length scales enter GL-equations on equal footing, 
however, any of normalizations might be considered as 
preferable, because the corresponding ^-dependencies in 
different normalizations look sometimes differently.) 

For small ip 2 / x 2 <C 1 one can drop the last term in 
Eq. (1) [which accounts for the screening current], so 
U = Uo(x) = —m + h^x 2 /2, B = H (no screening), and 
the system (1),(2) reduces to a single equation 



dx 2 x dx 
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We shall denote the solutions of the nonlinear equation 
(6) as tp(x) = K m (x). The analytical form of the (quasi- 
Kummcr) functions K m (x) is not known, but they may 
be easily found numerically from Eq. (6) [see Sec. IV]. 

In the limit ip — > 0, one can drop the term ■0 3 from Eq. 
(6) and consider the linear equation [5] 
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with the solutions denoted as tp(x) = K m (x). The func- 
tion K m (x) may be written in the analitical form [10] 
K m (x) = y m / 2 e- y / 2 T{y) (where y = ~/x 2 , 7 = h 6 /2). 
The function F{y) satisfies the confluent hypergeometric 
equation 

yT" + y)T' - vT = (8) 

[where fi — m + 1 is a positive integer, v = (1 — h^ 1 )/2], 

with a general solution T = C\F + C2F. The function 
F and can be written as an infinite series expansion 
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[We use the notation F(v,/j,,y) [11-13] for the Kummer 
functions, instead of M(—n,l + 1,2/), used in [14,15]. If 
v is a negative integer, or zero, F reduces to the polinom 
of power v. The function F has a logarithmic singular- 
ity at x — [16] and sometimes is dropped (C2 = 0). 
However, the function T{x) has a pre-factor y m / 2 ~ x m 
(see above), which cancels this singularity, so in a general 
case C 2 7^ (for m > 0).] 

The second of the boundary conditions (5) is equivalent 

to 



= 0, 



(9) 



x—R^ 



which is a nonlinear equation for v = (1 — h^ 1 )/2 (if m 
and i?£ are fixed). From Eq. (9) the proper value of 
fi£ (and also the space-profile ip(x)) may be found. In a 
general case, Eq. (9) has two roots for fi£ (if m > 0), what 
corresponds to two different solutions for ip(x) (see Sec. 



IV). The maximum field h^, above which there exists only 
the solution ^(x) = 0, corresponds to the upper critical 
field h 2 . 

[The solutions, found from Eqs. (l)-(9), describe the 
radially symmetric giant-vortex states of fixed vorticity 
m. The more complicated multi-vortex solutions (with 
the same total vorticity m = Sjmj, where nii is the vor- 
ticity of a vortex, situated at some point on the cylin- 
der cross-section) should be described by the equations 
in partial derivatives and will not be considered in the 
present paper. They may be studied numerically by the 
methods, used, for instance, in [17-19] to describe the 
results of experiments [20] with thin mesoscopic discs.] 



III. THE PHASE VIEWS 

The solutions of Eqs. (l)-(5) depend on the space co- 
ordinate x and several parameters [for instance, tp( x ) = 
i/)(ni, R,£, x, x); analogously for the potential U(x) 
and the field b(x)\. In this Section a general view is given 
of the solutions structure in this many-parameter nonlin- 
ear problem. Some details of this general picture will be 
discussed in the ensuing Sections. 



A. Field increase regime 

Let the vorticity m be fixed (m = 0, 1, 2, . . .) and con- 
sider at first the case m — (the vortex-free Meissner 
state). Consider the plane of parameters(i?£, x) (Fig. 
1(a)). To every point of this plane corresponds a set of 
solutions of Eqs.(l)-(5), which depend paramctrically on 
the external field (or h\ = x 2 h^). This set of solutions 
is unique for each point (R^ , x) and may be character- 
ized, for instance, by the field dependence of the maxi- 
mum value of the order parameter ip m ax(h^) in this point. 
It is helpful to look at the plane (R^, x) from above, and 
to imagine a peep-hole, pierced in arbitrary point of this 
plane, which allows to see the corresponding dependence 
V'mox (Such peep-holes are depicted as small circles 
in Fig. 1(a).) In doing so, one can find three regions in 
Fig. 1(a), marked as sf\ and Sjj. [The upper index 
A denotes a possibility of the hysteresis in these regions 
(see below). We shall mention Fig. 1(a), depicted on 
the plane (i?£, x) (and similar ones, like those in Fig. 3), 
as the phase views, to distinguish them from the phase 
diagrams, depicted on the plane (R^, h^) (as in Fig. 5).] 

Fig. 1(b) shows (schematically) the behavior ipmaxh^) 
in the region (the solid curve, the field ti£ increases 
from = 0). For small h^, the field is almost com- 
pletely expelled from the superconductor (the Meissner 
state, m = 0, 6 = 0). This part of the solid curve does 
not depend on and may be described by the London 
theory [2] approximation (tjj = 1). In larger fields, the 
order parameter depends on rather weakly. This part 
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of the curve may be described by the Ginzburg approx- 
imation [3] (ip = VK^i) = const). In still larger fields, 
a full system of GL-equations is needed to describe the 
behavior ip(h^;x). There exists a critical field (hi), at 
which the stable Meissner-like solution (the upper sec- 
tion of the solid curve) becomes absolutely unstable rel- 
ative small perturbations in its shape, and the solution 
ip(h^;x) acquires new stable space-profile, which may be 
called as the edge-suppressed state [7] (see the lower part 
of the solid curve, which is labeled also as "tail"; the 
difference between the various forms of solutions is illus- 
trated in Fig. 2). The transition from the upper to lower 
solid branch in Fig. 1(b) is a first-order phase transition 
(with a jump Si). For fields > hi the order parameter 
diminishes gradually and vanishes finally at the field hi 
by a second-order phase transition. Evidently, for fields 
/i£ <~ hi (when ip max <C 1) the linear equation (7) is 
applicable. 

The width of the tail (w = hi — hi), where the su- 
perconducting state (m = 0) is destroyed by the second- 
order phase transition to normal state, diminishes with 
the diminishing radius R^. The critical radius, at which 
the tail vanishes (w = 0), depends on x, and is depicted 
in Fig. 1(a) by a solid line. This line (<Si_n) divides the 
regions (s^) of first and (sjj) second-order phase tran- 
sitions. In the region sf" (marked as "no tails"), the 
transition to normal state is always of first order, by a 
jump from a finite value ip m ax to ip = 0. 

Type-I superconductors are naturally defined in our 
case as those, which belong to the region sf\ and type-II 
- to the regions s^, or Sjj. We see, that the boundary 
between type-I and type-II superconductors is a compli- 
cated function of (R^ , x) , it depends on the specimen 
geometry, and differs [7] from a simple phase boundary 
x = l/\/2 [8], which describes the case of an infinite 
(open) superconducting system. For instance, the super- 
conducting cylinder with x = 1 changes its type from 
type-II (above the point a in Fig. 1(a)) to type-I (be- 
tween the points a and (3), and again to type-II (below 
the point /?). 

B. Field decrease regime, type-II region 

If the field h^ decreases, starting from the large val- 
ues (h > hi, ip = 0), a weak superconducting state 
(with 'ipmax "C 1), which nucleates at hi, grows grad- 
ually into the edge-suppressed form (with ip max < 1 at 
/i£ = hi), but for h^ < hi the edge-suppressed solution 
transforms continuously into new "depressed" form (rep- 
resented by the dotted curve in Fig. 1(b)). In the field 
interval A = hi — h r there exist simultaneously two sta- 
ble solutions of Eqs. (l)-(5) (the solid and dotted lines 
in Fig. 1(b)), the depressed solution being responsible 
for the hysteresis in type-II cylinder. At the restoration 
field h r the depressed solution becomes absolutely un- 



stable and transforms by a first-order jump (5 r ) into the 
Meissner-type form (a solid line; the space-profiles of var- 
ious solutions are illustrated in Fig. 2). The width of the 
field interval, A = hi — h r , where hysteresis is possible, 
diminishes with the diminishing radius R^. The critical 
radius, at which A = 0, is presented by the dashed line 
in Fig. 1(a). No hysteresis is possible below this line. 

The first-order jumps (S r and Si in Fig. 1(b)) also di- 
minish with radius, they disappear simultaneously at the 
line A = (Fig. 1(a)). At this line the two disconnected 
pieces of the solid curve (Fig. 1(b)) met, and it acquires a 
continuous form, with an inflection point (see Fig. 1(d)). 
This curve is reversible, in the sense, that for each h^ 
there exists only one solution, which is independent of 
the field regime. For even smaller radius R^ (in the re- 
gion Sjj, below the line A = in Fig. 1(a)) the order 
parameter ip m ax(h^) behaves as is shown by the dashed 
line in Fig. 1(d). There is no hysteresis in the region Sjj, 
and the transition to normal state is always of second 
order, independently of x- value. 

C. Type-I region 

The order parameter ip m ax(h^) behaves very differently 
in the type-I region of Fig. 1(a). If the field increases, 
and if the radius R^ is sufficiently large, and the parame- 
ter x is sufficiently small, there are no tails in (see Fig. 
1(c)). The first-order jump (Si) to the normal state oc- 
curs at the point hi (where the superconducting solution 
in the field increase regime becomes absolutely unstable). 
For fields h^ > hi the normal state solution in Fig. 1(c) 
is absolutely stable, there are no other solutions here, but 
ip = 0. 

If the field h^ decreases now below hi, the supercooled 
normal state solution (ip = 0) remains stable relative 
small perturbations with ip ^ 0, down to the point h p . 
For fields inside the interval h p < h^ < hi there are two 
stable solutions for ip: one is ip = 0, the other (ip <~ 1) 
is represented by the upper section of the solid curve 
in Fig. 1(c). (Both solutions are stable, because small 
perturbations in there shape die down.) Thus, in the 
field interval A = hi — h p the supercooling of the normal 
state is possible [see the hysteresis loop in Fig. 1(c)]. 

Below h p the normal state becomes absolutely unstable 
(small perturbations grow), and a feeble superconduct- 
ing solution (or, the "precursor" state, ip -C 1) establishes 
in the bulk of a cylinder. (The position of h p coincides 
exactly with the critical field of x-independent second- 
order (s, n)-transitions and may be found from a linear 
theory [5].) The precursor solution exists down to some 
restoration field h r , where it becomes absolutely unsta- 
ble and the Meissner form restores [notice the presence 
of first-order jump S r in Fig. 1(c)]. The field interval 
A p = hp — h r , where the precursor solution exists, is 
usually very small (A p ~ 10~ 2 — 10~ 4 h p ), so in type- 
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1 superconductors the restoration of superconductivity 
from a supercooled n-state is "almost first-order" phase 
transition. The exact value of the restoration field (h r ) 
and the corresponding amplitude of the precursor solu- 
tion (ip r ) can not be found from the linear theory, but 
self-consistent solution of the full system of nonlinear GL- 
equations is required. 

If the radius R^ diminishes, the jumps S r and Si also 
diminish (as well as the interval A = hi — h p , see the 
dashed curves in Fig. 1(c)), they all disappear simulta- 
neously with A — > 0. The line A = (for small x and 
merges with the line Si_n in Fig. 1(a) at the point 
(3 (compare the dashed curves in Figs. 1(c) and 1(d)). 
No supercooling of normal state is possible below the line 
A = 0. 

D. Examples of co-ordinate dependencies 

The co-ordinate dependence of the solutions [as they 
are seen through the peep-hole, pierced in the type-II 
point i?f = 10, >f = 1.2 in Fig. 1(a)] is shown in Fig. 2 
for different values of the external field H (h^ = H/H^). 
In the field increasing regime, the sequence of curves, 
which result from Eqs. (l)-(5), are numerated as (0, 
1, 2, 2e, K m ). The order parameter ip(r) (Fig. 2(a)) 
changes its form continuously from (tp = 1 at h^ = 
0) to 1 (h s = 0.9037), and then to 2 {h t = 1.0951 - 
this is the critical field, marked as hi in Fig. 1(b)). At 
this point the solution 2 becomes absolutely unstable and 
(by a first-order jump, <5i = 0.155) acquires the edge- 
suppressed form 2e (h^ = 1.0952). If the field increases 
further, the solution 2e passes continuously into the form 
K m (the Kummer-type solution), which vanishes finally 
at the critical field h^ = h 2 = 1.0016, where ip = 0. 
Fig. 2(b) illustrates the behavior of the induction B{r) 
versus the external magnetic field H. At the point hi the 
solution 2 switches (by a jump <5i) from the Meissner form 

2 (the external field is effectively screened out) to the 
edge-suppressed form 2e (the external field penetrates the 
edge region without screening). The curve K m may be 
described by the linear equation approximation (B = H) . 

The sequence of solutions, appearing in Fig. 2 in 
the field decreasing regime, is different (see Fig. 1(b)). 
The normal state solution (ip = 0) is absolutely sta- 
ble for /i£ > h 2 , and absolutely unstable for < h 2 . 
If the field decreases, a small superconducting solution 
(of the Kummer-type, K m ) appears at — h 2 and 
transforms continuously into the edge-suppressed form 
2e (with ip ~ 1 at = h 2 ). If the field decreases be- 
low hi, the solution 2e keeps transforming continuously 
into the form Id (at h^ = h r ). The intermediate sequence 
of curves (2e <==^> Id) presents the depressed solutions, 
which arc characterized by a smaller average values of the 
order parameter ip, in comparison with the Meissner so- 
lutions (1 ^=^> 2). For h^ < h r the depressed solution Id 



becomes absolutely unstable and transforms (by a jump 
S r at h^ — h r ) into the Meissner form 1. It is important, 
that in the interval of fields A~hi — h r there exist simul- 
taneously (for the same h^) two independent solutions of 
Eqs. (1) (5) (the Meissner and depressed forms). The 
presence of two solutions means a possibility of hysteresis 
in the system. 

[The Meissner solution can be obtained from Eqs. (1) 
(5) by the iteration procedure [21], started with a large 
trial function tp ~ 1. The depressed form is obtained by 
the iterations, started with a small trial function ip ~ 
0.01. Outside the interval A = hi — h r both iteration 
procedures produce the same self-consistent result.] 

Additional examples of co-ordinate dependencies of the 
solutions for various R and x may be found in [7] . (The 
space-profiles of the precursor solutions are depicted be- 
low in Fig. 8.) 

E. The phase views for m=l,2 

The analogous phase views on the plane (R^,x) for 
the vortex states m = 1 and m — 2 are depicted in Figs. 
3(a,b), they are similar to the phase view m = in Fig. 
1(a). The only essential difference is in the presence of 
the minimal radius R^ [see the dashed lines C sn in Figs. 
3(a,b)], below which only the normal state solutions are 
possible. This is natural, because in a case of very small 
radius specimens the intrinsic magnetic field of a vortex is 
too strong for the survival of the superconducting state. 

Fig. 4 illustrates the solutions behavior in the vortex 
states m = 1, 2, as they are seen through the peep-holes 
in type-II points R^ = 10, x = 1.2 in Figs. 2(a,b). If the 
field increases, the sequence of the soluions is (0, 2, 2e, 
K m ). In the Meissner-type states (0, 1, 2) the external 
field is screened out. At the field h^ = hi the jump trans- 
formation to the edge-suppressed form (2-^2e) occurs 
(with a jump 5i). (The concrete values h r ,hi,h 2 ,S r , Si 
are given in the caption to Fig. 4.) 

In the field decreasing regime the sequence of the ap- 
pearing solutions is (K m , 2e, Id, 1, 0), with a jump 
transformation (S r ) from Id to 1. In the field interval 
A = hi — h r there are two solutions: one is the Meissner- 
type form (intermediate between 1 and 2) , and the other 
is the depressed form (intermediate between 2e and Id). 
Within the field interval A the hysteresis transitions be- 
tween the Meissner-type and depressed states of the same 
vorticity m are possible. For fields outside the interval 
A no hysteresis is possible, because only one solution 
exist there (see Fig. 1(b)): either the Meissner-type so- 
lution (for h^ < h r ), or the edge-suppressed solution (for 
h ( > hi). 

The giant vortex states with m > 2 may be studied 
analogously. 
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IV. THE PHASE DIAGRAMS 

In this Section the comparison is made of the phase 
diagrams (or, the critical fields h^) for the second-order 
(s, n)-transitions in the increasing magnetic field, found 
sclf-consistently and in the linear approximation. [See 
also Ref. [7a], where the critical fields h\ for the first- 
order jumps to the edge-suppressed states were studied.] 

From the phase view of Fig. 1(a) it follows, that if the 
magnetic field increases, the order parameter vanishes 
either gradually [with a tail in the curves of Figs. l(b,d)], 
or by the first-order jump [Fig. 1(c)]. Fig. 5(a) is a phase 
diagram, which shows the dependence of the critical field 
h c sn in Fig. 1 [i.e. the field h\ (type-I), or hi (type-II), 
for which the transition from s- to n-state occurs (in the 
field increase regime)], as a function of R^, for different 
x and to = 0. The thick line K m corresponds to the 
(s, n)-diagram, found from the linear equation (7) with 
the boundary conditions (5). 

As is seen from Fig. 5(a), for x > 1 the self-consistent 
critical curves h c sn {h{) coincide with the universal curve 
K m {h^) (found from the linear equation). Thus, for type- 
II superconductors (with x > 1) the linear approach [4,5] 
gives correct description of the (s, n)-boundary. 

However, for type-I superconductors (with x < 1), the 
critical curves h c sn {h^) deviate from the universal curve 
K m (h^) considerably. The inspection of the asymptotes 
of the curves h° sn (for ^> 1) shows, that they have dif- 
ferent functional dependences versus parameter x. For 
type-II superconductors (x > 1) the asymptotes are 
h c sn « 1 (they are x-independent, in accordance with the 
linear theory [5]). For type-I superconductors (x < 1) 
the asymptotes behave as h c sn « x~ p (1 < p < 2); this 
means, that the universal linear theory fails for small 
x < 1 and large R^ . (For sufficiently small radius R^ the 
(s, n)-transition is always of second order, so in this case 
the linear theory is valid for all x.) 

The analogous conclusions may be drawn from Fig. 
5(b) (to = 1) and Fig. 5(c) (to = 2). 

The fact, that the critical fields h c sn in Fig. 5 for x < 1 
deviate from the predictions of the linear theory [5] , can 
be attributed to strong nonlinear competition between 
two characteristic lengths A and £ in GL-equations. For 
x = A/£ > 1 the nonlinear equations can be linearized, 
but for x = A/£ < 1 they can not. (Evidently, the lim- 
itations for the linear theory can not be deduced from 
the linear theory itself, but self-consistent analysis of full 
nonlinear system of equations is required.) 

A. Two solutions of the linear equation 

Notice the peculiar behavior of the phase diagrams in 
Figs. 5(b,c). For sufficiently small R^ the line R^ =const 
crosses the bottom part of the curves K m at two points. 
This means the existence of two solutions K m (x) of the 



linear equation (7) , which both satisfy the boundary con- 
ditions (5) for the same R^, but with different fi£ (see Eq. 
(9)). These two solutions are depicted in Fig. 6(a) (the 
solid lines 1 and 2), as they are seen through the peep- 
holes in Fig. 5(b), pierced along the line R^ = 1.5 at the 
points h,£ = 0.67 and h 6 = 2.32 (to = 1). 

[The insert to Fig. 6(a) illustrates (schematically) the 
behavior of the solutions of the hypergeometric equation 
(7) for to = 1. The function K m (x), in a general case, 
has two extremal points A and B, where K' m (x) = 
(marked by the arrows), whose positions depend on h^. 
The curves A and B in Fig. 6(a) have different functional 
behavior: the curve A grows monotonously, reaching the 
extremum at x = R^ = 1.5. The curve B has a zigzag in 
the vicinity of x = R^ (the zigzag amplitude S is small 
for small R^ and is not seen in this scale). The functions 
K m (x) for to > 1 behave analogously, see Fig. 6(b) for 
to = 2, R^ = 2.5.] 

The existence of two different solutions of Eq. (7) has 
clear physical interpretation. In the absence of the exter- 
nal field, the vortex is held inside the cylinder by pinning 
to the boundary (which is the source of the inhomogene- 
ity in otherwise homogeneous system). However, in small 
radius cylinders the internal vortex field Bi is too strong 
to be confined inside the mesoscopic sample by the pin- 
ning force, and it breaks outside. To prevent the field Bi 
from leaking to the outer space, it is necessary to impose 
a finite external field H, which helps to keep Bi inside 
the specimen. (This can be considered as an example of 
the so-called re-entrant superconductivity, or, the field- 
enhancement effect.) However, if the field H increases, 
the superconductivity will be finally destroyed. Two so- 
lutions (A and B in Figs. 6(a,b)) describe two different 
physical situations. The solution B corresponds to the 
superconducting vortex state (to > 0) being destroyed 
by a large external field H . The solution A (with smaller 
H) corresponds to the vortex state being destroyed by 
the internal field of it's own vortex. 

The solutions A and B in Figs. 6(a,b) are depicted 
as they are seen through the peep-holes in Figs. 5(b,c), 
which lie near the opposite sides of the thick curves K m 
[these curves represent the states with ip max *C 1; the cor- 
responding solutions are the hypergeometric func- 
tions K m (x) of the linear Eq. (7)]. In the intermediate 
peep-holes i (Fig. 5) the solutions ip(x) have finite am- 
plitude ipmax, but (if tpmax is yet sufficiently small) the 
solution may be approximated by the function K m (x) of 
the nonlinear Eq. (6) [see the dashed curves C in Figs. 
6(a,b)] . To find ip(x) with still larger ip ma x, it is necessary 
to solve the full system of Eqs. (1),(2). Such solution is 
presented by the curve D in Fig. 6(b). 
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B. Supercooling of the normal state 

As was explained in Sec. Ill, in any peep-hole in the re- 
gion of Fig. 1(a) (wherex > 1 and A > 0) there exist 
two solutions of Eqs. (1),(2): one is the Meissner-type 
solution, and the other is the depressed (or partially sup- 
pressed) solution, which is responsible for the hysteresis 
in type-II superconductors (see Fig. 2 for m = and Fig. 
4 for m = 1, 2). 

The analogous hysteresis phenomena exist also in type- 
I superconductors [in the region s^ of Fig. 1(a), where 
K < 1 and A > 0; analogously for Fig. 3]. In this 
case, there are also two solutions: one is again of the 
Meissner type, but the partially depressed solution (rep- 
resented by the dotted curve in Fig. 1(b)) degenerates 
now into the totally depressed normal state of Fig. 1(c). 
Thus, in type-I superconductors the second branch of so- 
lutions in Fig. 1(b) corresponds to the supercooled n-state 
(ip(x) = 0). (Notice, that in typc-II superconductors the 
supercooled normal state can not exist.) 

The curve sh in Fig. 7(a) corresponds to a maximal 
field, at which the superheated Meissner state (m = 0) 
may still exist in type-I superconductor with x = 0.1 in 
the field increase regime (i.e. the field ft, 2 in Fig. 1(c)). 
The curve sc corresponds to a minimal field, at which 
the supercooled normal state may still exist in the field 
decrease regime (i.e. the field h r in Fig. 1(c)). These two 
curves merge at the point LCP (the Landau critical point 
[22]) into a single curve (the calculated difference between 
sc- and sft-curves at LCP is less than 1 • 10~ 4 ). The free 
energies of the superconducting states, which lie between 
sc- and sft-curves, depend on the field h%, and at some 
value of fi£ the difference of free energies AG = G s — G n 
vanishes. The curve eq in Fig. 7(a) shows (schematically) 
the position of the equilibrium curve (AG = 0). All three 
curves merge at LCP, which is a three-critical point. 

The analogous dependences are depicted in Fig. 7(b) 
for x = 0.5. It is clear, that the width of the region 
between sc- and sft-curves diminishes with x increasing; 
for x = 1 both curves merge into a single curve and the 
(normal-state) hysteresis region disappears. 

Notice, that for R^ 3> 1 the maximal supercooling 
of the normal state in type-I superconducting cylin- 
der is reached in the field fi£ — 1, i.e. in the field 
= H C 2 = <fro / (2tt£ 2 ) . To the same conclusion came 
Ginzburg [3] and Abrikosov [8], who used the approxi- 
mations ifr = const and [The dependence of the 
superheating field H s h (H) for type-I cylinders was found 
earlier (in a different paramctrization) by Esfandiary and 
Fink [23] , who used self-consistent solutions of the nonlin- 
ear equations (the supercooling field H sc was not studied 
in [23]). (See [6] for the review of early theoretical and 
experimental works on the problem of hysteresis in type-I 
superconductors.)] 

The comparison of sc-curves in Figs. 7(a,b) with the 
universal K m -cmve in Fig. 5(a) reveals, that they are 



identical and have the following asymptotic behavior: 
h t (Re) = 1 for i? c > 1; Rt(h s ) = 2.8/ h t for ft € > 2 
(see dotted lines a in Figs. 7(a,b)). Notice, that these 
curves coincide with each other for all and x < 1 (not 
only for R^ ^> 1 and x <C 1, as was found in [3,8]). This 
coincidence is not accidental, because at the supercooling 
n-boundary the precursor solutions appear (ip max *C 1) 
which are described by the x-independent linear equa- 
tion (7). 

The precursor solutions are depicted in Fig. 8 for the 
fields fi£, which just precede the jumps to the Meissner 
state with ip ~ 1. The precursor solution shape is x- 
independent (in accordance with Eq. (7)), but the ex- 
act position of the jump and the amplitude of ifr ma x are 
x-dependent and should be found from the full system 
of GL-equations (1), (2). Notice, that the precursor so- 
lutions describe the superconductivity nucleation in su- 
percooled type-I cylinder as a bulk (though very feeble) 
effect. [The detailed study of the precursor solutions be- 
havior would be presented elsewhere.] 

C. The free-energy functional 

As was mentioned above, the problem of supercool- 
ing and superheating in type-I superconductors was con- 
sidered first by Ginzburg [3], who analysed the behav- 
ior of the free energy functional in the approximation 
ifr = const (m = 0). It is expedient to compare the re- 
sults of the self-consistent and approximate approaches. 

In Fig. 9(a) the exact dependence ip m ax(h\) is shown 
for x = 0.5 at R^ = 3 and R^ = 1 (m = 0). It is evident, 
that the width A of the hysteresis region diminishes with 
i?£ diminishing, and vanishes at some R m in (at the point 
LCP in Fig. 7(b)). 

In Figs. 9(c) the exact field dependence is shown of 
the normalized free energy (Ag) for x = 0.5 and R^ = 3. 
(The expressions for Ag and magnetization (— AttM) may 
be found in [7]). Evidently, the part of the curve Ag(h^), 
which lies to the right of the point eq (where Ag = 0) , 
corresponds to the metastable s-state with Ag > 0. In 
this region s-state is energetically unstable (because n- 
state has smaller free energy, Ag = 0). [However, the 
energetically unstable s-state is stable relative small spa- 
tial perturbations.] To the left of the point eq the s-state 
is energetically more favorable (Ag < 0) , and the n-state 
(with Ag = 0) is metastable. [However, the energetically 
metastable n-state is stable relative small spatial distur- 
bances.] At the equilibrium point eq the transition from 
one branch of the solution to the other (with smaller free 
energy) may happen. In this case the system behavior 
would be totally reversible (without hysteresis). 

In Fig. 9(b) the field dependence of the magnetization 
(— 4nM) is illustrated. In the case of equilibrium tran- 
sition the reversible jump of the magnetization should 
be observed at the point eq. If the metastable states 
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are realized, the jumps of the magnetization may hap- 
pen anywhere between the points sc and sh, with the 
accompanying hysteresis. [There is a principal difference 
between the jump transitions in Fig. f (b) and Fig. 8. 
The first-order jumps in type-II superconductors (Fig. 
1(b)) occur within the same state (to = 0), while the 
jump transitions in type-I superconductors (as in Fig. 9) 
occur between different states (s and n).] 

Notice, that in the approximation tp = const [3] the 
free energy AG(ip) has the form, schematically depicted 
in Fig. 10(a) (the curves I-t-5 correspond to the increas- 
ing fields H). The extremas of these curves define the 
possible values of the order parameter tpo(H), which are 
depicted (schematically) in Fig. 10(b). [The numerals 
l-i-5 in this figure mark the positions of the extremas 
on the corresponding curves in Fig. 10(a).] Thus, in 
the approximate approach [3] there arc two solutions for 
ipo 7^ in the field interval H sc < H < H s h, one (which 
corresponds to the minimum of the free energy AG) is 
energetically stable, the other (which corresponds to the 
maximum of the free energy) is energetically unstable. 

However, according to the self-consistent theory, in 
type-I superconductors there exist only one solution 
with ip ^ [which might be energetically stable, or 
metastable, depending on H (see Fig. 9(a))], the other 
(with tp = 0) corresponds to the supercooled normal 
state. Thus, according to the exact theory, the unsta- 
ble branch ip ^ of approximate theory does not exist 
at all. 

This contradiction arises because of the following rea- 
son. In the rigorous theory the free-energy is a func- 
tional, G[ip{x)\, which produces a number, G, from a 
function, ip(x). The function ip(x) depends parametri- 
cally on H, so the functional G[ip(x)] is a function of H. 
In the Ginzburg approximation the functional G[ip(x)] 
is replaced by a function G(ipo,H), where ipo is consid- 
ered as an arbitrary number. This function has minima 
and maxima, what is illustrated by Fig. 10. However, 
as follows from the self-consistent solution of nonlinear 
problem, some of the values ipo are forbidden. As can be 
seen from Fig. 9(c), the exact functional G(H) increases 
with H, but it does not reach an extremal point, where 
G+(H) = 0, and drops to zero at some value of H, where 
ip(H) terminates. Thus, the Ginzburg approximation for 
G(V>o) [3] is qualitatively valid, if tp(x) w ipo = const 
[see the stable (solid) branch of the curve ipo(H) in Fig. 
10(b)], but fails for those (forbidden) values of ipo, which 
belong to the unstable branch of ipo (the dashed curve). 
Nevertheless, the point of maximal supercooling, H SCl 
is described by the Ginzburg approximation correctly, 
because at this point (ipo = 0) the precursor state nu- 
cleates. [The precursor solution has very small ampli- 
tude, ip(x) ~ 0: what always satisfies the condition [3] 
ip(x) w const, independently of the real space-profile of 
ip(x) (see Fig. 8).] 

We shall restrict ourselves by these methodical re- 



marks, while comparing the results of the rigorous and 
approximate approaches. 

V. CONCLUSIONS AND DISCUSSION 

We have studied in detail the one-dimensional solu- 
tions of GL-equations in a case of cylinder geometry. The 
main new results of the present investigation are sum- 
marized below. Mention among them: the phase views, 
presented in Figs. 1(a) and 3(a,b); the existence of the 
edge-suppressed and depressed solutions (Figs. 1,2,4); 
the phase diagrams of Fig. 5; the discussion of the ap- 
plicability of the linear equation approximation (Sec. V); 
the discussion of two independent solutions of the lin- 
ear equation (Fig. 6); somewhat novel insight into the 
problem of hysteresis in type-I superconductors (Figs. 7- 
9) and the discussion of the previously unknown precur- 
sor solutions; the discussion of the free-energy functional 
and the comparison of the rigorous and approximate ap- 
proaches to the hysteresis problem (Fig. 10). 

In conclusion, some additional comments to the topics, 
discussed above, should be made. 

The edge-suppressed solutions in a cylinder geometry 
were described first by Fink and Presson in a footnote to 
their paper [9], but were disregarded as being unstable 
in the energy space. We consider these solutions (as well 
as depressed and precursor solutions) as stable in the co- 
ordinate space. Mathematically, they are usual axially 
symmetric solutions, which have different forms in dif- 
ferent regions of parameters (due to the nonlincarity of 
the problem). These solutions describe symmetrical in- 
states, which the physical system may occupy. At the 
critical field hi the Meissner state (to = 0) becomes un- 
stable and the external field fluds the outer region of a 
cylinder (see the edge-suppressed solution 2e in Fig. 2). 
Another possible mechanism for the field penetration is 
a creation of a vortex line (to = 1) on the cylinder axis 
(see solutions in Fig. 4). Comparing the free-energies 
of the competing states, one can find the points of equi- 
librium transitions [7], and also the field interval, where 
the edge-suppressed state may exist as a metastable one. 
The metastable states may manifest themselves in exper- 
iments, where such metastability is realized (for instance, 
in the hysteresis phenomena [24] , in paramagnetic Meiss- 
ner effect [25], etc.), especially in a case of mesoscopic 
samples (with few number of vortices), when the stabi- 
lizing role of the boundary is important. 

It is interesting, that similar effects - such as jumps of 
the magnetization within the states of fixed to, the tran- 
sition to the edge-suppressed form of a giant- vortex solu- 
tion, the re-entrance of superconductivity and the field- 
enhancement effects - exist also in a case of mesoscopic 
disks of small height, studied self-consistently in [17,18]. 
This means, that these nonlinear effects are common to 
various sample geometries and are determined mainly by 
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the parameters, entering the GL-equations, so, the role 
of the boundary is, probably, not crucial. 

No immediate comparison between the present the- 
ory and experiment was attempted, partly because the 
model case of infinitely long cylinder approximates only 
remotely the thin-disk geometry mostly used in recent 
experiments [20]. Moreover, in discussing such experi- 
ments, it is necessary to consider also the asymmetric 
multi- vortex solutions (see the end of Sec. II), what re- 
quires using specific numerical methods and large com- 
puters (see, for instance, [17,18]). However, a num- 
ber of qualitative predictions, obtained from the one- 
dimensional theory, may be used to interpret some of 
the peculiarities, observed on mesoscopic samples. For 
instance, it might be possible to attribute some jumps 
of the magnetization (seen in [20,24,25]) not to transi- 
tions between different m-states, but occurring within 
the same m-state (due to the order-parameter reconstruc- 
tion, while it passes to the edge-suppressed form). How- 
ever, a special experimental analysis would be needed to 
confirm the existence of such transitions. 

Notice, that the precursor solutions, studied above (see 
Fig. 8), show, that type-I superconductor may be found 
in two different (metastable) hysteretic states (in addi- 
tion to the stable Meissner state). One is a supercooled 
normal metal (tp = 0); this state may persist down to 
the maximal supercooling field (h p ). The other is the 
supercooled precursor state (ip 0), which nucleates at 
hp and survives down to the restoration field (h r ), when 
the first-order transition to the Meissner state (ip » 1) 
occurs. This precursor state is more prominent for val- 
ues of x, laying in the vicinity of S , i_n-phase boundary 
(Fig. 1(a)). Though the field interval A p , where the pre- 
cursor state may be found, is small, it would also be of 
interest to seek for experimental evidences of its possible 
existence. 

Evidently, further theoretical and experimental study 
of the problems, discussed above, is desirable. 

[The one-dimensional solutions in other geometries (a 
plate of finite thickness and superconducting half-space) 
were addressed earlier in a number of publications (see, 
for instance, [26-30]). This topic will be discussed else- 
where.! 
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Figures captions 

Fig. 1. (a) - The phase view for a cylinder in the 
vortex-free Meissner state (m = 0) in the field increasing 
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regime; Sj - the region of the first-order phase tran- 
sitions from s- to n-state; and Sjj - the regions of 
second-order phase transitions. The curve Si_n divides 
the regions of first- and second-order phase transitions. 
In the region above the dashed line A = the hysteresis 
transitions are possible; no hysteresis is possible below 
this line, (b) - Schematic behavior of the order param- 
eter ip m ax(h^) in the region (m = 0, A > 0, w > 0). 
(c) - Schematic behavior of the order parameter in the 
region (m = 0, A > 0, w = 0). (d) - Schematic be- 
havior of the order parameter in the region Sjj (to = 0, 
A = 0, w > 0). The used notations are explained in the 
text. 

Fig. 2. The coordinate dependences: (a) - for ip(x) 
and (b) - for b(x), as they are seen through the peep- 
hole in Fig. 1(a) at R 6 = 10, x = 1.2 (m = 0). The 
curves are calculated for the dimensionless fields = 
H/Hf. - h£ = 0; 1 - ht = 0.6275; Id - h 6 = 0.6276; 
2 - h(. = 0.7605; 2e - h ( = 0.7606; K m - h ( = 0.992; 
i'max = at /ij = 1.000. The jump Si of i\) max (see Fig. 
1(b)) at the transition ld-^1 is Si = 2 • 10~ 4 . The jump 
at the transition 2-^2e is S r = 0.155. The arrows show 
how the solutions transform, if the field is increased, 
or decreased. The depressed solutions, responsible for 
the hysteresis, have the form, intermediate between 2e 
and Id. 

Fig. 3. The phase views: (a) - for m = 1 and (b) - for 
to = 2. Notations are the same, as in Fig. 1(a). Below 
the line C sn lies the normal-metal region n. 

Fig. 4. The coordinate dependences i[)(x) and b(x), 
as they are seen through the peep-holes in Fig. 3(a) at 
Rt = 10, x = 1.2 (to = 1), and in Fig. 3(b) (to = 2). 
The curves in (a) and (b) are calculated for the fields 
hf. - h s = 0; 1 - h£ = 0.6367; ld-h^ = 0.6368; 2 - 
h ( = 0.7610; 2e-h i = 0.7611; K m - h 6 = 0.986. The 
jump at the transition ld-^1 is Si = 0.018; the jump 
at the restoration transition 2—>2e is A r = 0.238. The 
curves in (c) and (d) are calculated for the fields h^: 

h 6 = 0; 1 - h 6 = 0.6358; Id - h 6 = 0.6359; 2 - h t = 
0.7619; 2e-h i = 0.7620; K m - h 6 = 0.986. The jump 
at the transition ld—*l is Si = 0.030; the jump at the 
restoration transition 2-^>2e is S r = 0.244. 

Fig. 5. The phase diagrams represent the critical fields 
Ii£ = H/H^, at which the transition from s- to n-state oc- 
curs for a given = R/£ and various x (see the numer- 
als at the curves) and different vorticities to: (a) - to = 0; 
(b) - to = 1; (c) — to = 2. Thick curves K m correspond 
to the solutions of linear equation (7) . The curves with 
x > 1 are well represented by _ftT m -curve [thus, the linear 
theory [5] gives correct description of (s, n)-boundary] . 
However, for x < 1 (and large R) the curves deviate 
strongly from the curve K m [thus, for x < 1 the lin- 
ear theory [5] is not applicable]. For sufficiently small 
i?£ (when (s, n)-transition is of second order) the linear 
theory is valid for all x. 



Fig. 6. Two solutions (1 and 2 ) of the linear equation 
(7) in the case: (a) -to = 1 [R^ = 1.5, the fields fi£ = 0.67 
(1 ) and h/: = 2.32 (2 )]; (b) -m = 2[R i = 2.5, the fields 
h(_ = 0.291 {1 ) and = 1.742 (2 )]. The insert to 
Fig. 6(a) shows the schematic behavior of two possible 
solutions. The zigzag amplitude 5 for the solution 2 is 
not seen in this scale in (a), but is present in (b). (The 
arrows at the curves 2 show the maximums of ip(x).) 
The dashed curves 3 represent solutions of the nonlinear 
-ftT m -equation (6) in the intermediate peep-hole i [with 
fi£ = 1.5 (a) and ti£ = 1.6 (b)]. The dotted curve 4 
represents the self-consistent solution, seen through the 
peep-hole i' (hf. = 1) in Fig. 5(c). [All the curves in Fig. 
6 are normalized to ip max = 1.] 

Fig. 7. The upper critical field (the curve sh) for 
superheated s-state, and the lower critical field (the curve 
sc) for supercooled n-state of the vortex-free Meissncr 
state (to = 0) for x = 0.1 (a) and x = 0.5 (b). For 
x = 1 the sc- and s/i-curves practically coinside. The 
dotted curve a ~ 2.8/ gives good approximation to 
sc-curve for R^ < 2. The broken curve eq (schematic) 
corresponds to the equilibrium transition between s- and 
n-states. LCP is the tri-critical Landau point. The curve 
sc coincides with the phase boundary K m (dashed line), 
found from the x- independent linear equation (7). 

Fig. 8. Precursor solutions for R^ = 3.6, x = 0.9: (a) 
- ip(x), (b) - b(x). Supercooled precursor state nucleates 
at h p — 1.0231 (with ip w 0), and takes forms: 1 - at 
h 6 = 1.0222; 2 - at h s = 1.0210; 3 - at = 1.0123; 4 - 
at /i| = 0.9991 (this is maximally supercooled precursor 
state, which at h r = 0.9990 passes into the Meissner 
form). The dotted line l n is the normalized (Kummer) 
function J « 1), which can not be described as tj){x) ~ 
const [3]. However, the exact solution 1 is small (ip(x) w 
= const), so the field of maximal supercooling of n-state 
(hp) may be found from the Ginzburg [3] approximation 
(and also from the linear theory [5]). 

Fig. 9. The field dependence of the order parameter 
ipmax (a), the magnetization (M% = M/H£) (b) and the 
free energy Ag (c) in the Meissner state (to = 0) of the 
cylinder with x = 0.5 and R^ = 3. The curve i?£ = 1 in 
(a) demonstrates, that the hysteresis region A vanishes 
for small R^ < 1. 

Fig. 10. (a) - Free energy functional AG versus 
tp = const, according to the Ginzburg approach [3] 
(schematic). The sequence of curves 1-5 corresponds 
to the field H increasing. Minimums of AG correspond 
to stable states, maximums - to unstable states, (b) - 
The order parameter vp , found from the extremums of 
AG. Stable branch is depicted by solid line, unstable - 
by dashed line. The numerals 1—5 correspond to those 
in (a). The unstable branch does not exist in the self- 
consistent theory (see Fig. 9(a)). 
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